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ABSTRACT 

The design and later use of modern spectropolarimeters and magnetographs require a number of tolerance 
specifications that allow the developers to build the instrument and then the scientists to interpret the data 
accuracy. Such specifications depend both on device-specific features and on the physical assumptions un- 
derlying the particular measurement technique. Here we discuss general properties of every magnetograph, 
as the detectability thresholds for the vector magnetic field and the line-of-sight velocity, as well as specific 
properties of a given type of instrument, namely that based on a pair of nematic liquid crystal variable retarders 
and a Fabry-Perot etalon (or several) for carrying out the light polarization modulation and spectral analysis, 
respectively. We derive formulae that give the detection thresholds in terms of the signal-to-noise ratio of the 
observations and the polarimetric efficiencies of the instrument. Relationships are also established between 
inaccuracies in the solar physical quantities and instabilities in the instrument parameters. Such relationships 
allow, for example, to translate scientific requirements for the velocity or the magnetic field into requirements 
for temperature or voltage stability. We also demonstrate that this type of magnetograph can theoretically reach 
the optimum polarimetric efficiencies of an ideal polarimeter, regardless of the optics in between the modulator 
and the analyzer. Such optics induces changes in the instrument parameters that are calculated too. 
Subject headings: Sun: magnetic fields - Sun: photosphere - Sun: polarimetry 



1. INTRODUCTION 

Spectropolarimetry and magnetography have become two 
of the most useful tools in solar physics because they provide 
the deepest analysis one can make of light. Solar information 
is encoded in the spectrum of the Stokes parameters. We mea- 
sure this spectrum and infer solar quantities from it. Recently, 
less and less conceptual differences exist between spectropo- 
larimeters and magnetographs except for the specific devices: 
the formers usually include those instruments using a scan- 
ning spectrograph; the latter usually employ a bidimensional 
filtergraph like a Fabry-Perot etalon. Some decades ago 
magnetographs were only able to sample one or two wave- 
lengths across a spectral line; nowadays, new technologies 
provide a better wavelength sampling, thus enabling the sci- 
entists to interpret the data in terms of sophisticated inver- 
sion techniques of the radiative transfer equation, a procedure 
similar to the one regularly used with spectropolarimeters. 
Some of the instruments mentioned below enter this category. 
Modern solar spectropolarimeters and magnetographs are of- 
ten vectorial because all four Stokes parameters of the light 
spectrum are measured. Longitudinal magnetography (i.e., 
Stokes I ±V) can be interesting for some specific applica- 
tions, but the partial analysis is usually included (if possible) 
as a particular case of the more general, full-Stokes polarime- 
try. Some of these modern instruments have been recently 
built or are curre ntly in operation (e.g. , the T e nerife Infrared 
Polar imeter, TIP, IMartmez Pillet et al.l 119991 ICollados et all 
20071: the Diffraction-Limited S pectro-Polarimeter, DLSP, 



Sankarasubramanian et al.l 120031: the air-spaced F abry-Perot 
based CRISP inst rument. IScharmer et al.f [2008; the spec- 
tropolarimeter. SP.Lites et al.l 120011 for the Hinod e mission, 

~l 12001 the Visi- 
120 10t the Imag- 



iKosugi et alJ 120071 CRISP. iNaravan et" 



ble Imaging Polarimeter, VIP, iBeck et al 



ing M agnetograph experime nt, IMaX, Martinez Pillet et all 
HOilforthe Sun rise mission, iBarthol et al.ll201 ll and the He- 
lioseismic and Magnetic Imager, HMI, iGrahametafl I2003L 
for the Solar Dynamics Observatory mission, iTitld 120001) . 
Some other are being designed and built for near future op- 
eration and missions (e.g., the Polarim etric and Helioseismic 
Imager, SO/PHI, [formerly ca lled VIM, IMartmez Piliell2007H 
for the Solar Orbiter mission, iMarsch et al.ll2005l) . 

The common interest of users of these instruments is cen- 
tered in vector magnetic fields (of components B, 7, and </>) 
and line-of-sight (LOS) velocities (vlos)- Some spectropo- 
larimeters can provide information on temperatures as well 
(and eventually on another thermodynamical quantity) but 
that feature is not common to all of them. Therefore, assess- 
ing the magnetograph capabilities in terms of their accuracy 
for retrieving magnetographic and tachographic quantities is 
in order since such an analysis can diagnose how far reach- 
ing is our current knowledge of the solar dynamics and mag- 
netism. The diagnostics is relevant both for the design of new 
instruments in order to maximize their performances and for 
the analysis of uncertainties in data coming from currently 
operating devices. General considerations can obviously not 
be made but a few. Specifically, we here study in Sect. |2]the 
detection thresholds induced by random noise on the inferred 
longitudinal and transverse components of the magnetic field; 
in the particular case of photon-induced noise we also find un- 
certainty formulas. Both thresholds and relative uncertainties 
are obtained in terms of the signal-to-noise ratio of the obser- 
vations and of the polarimetric efficiencies of the instrument. 
Since such efficiencies vary from instrument to instrument, 
at that point, the analysis concentrates in a particular type of 
magnetograph, namely that consisting of two nematic liquid 
crystal variable retarders (LCVRs) as the polarization modu- 
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lator and a Fabry-Perot etalon. In Sect. [3] we demonstrate that 
these polarimeters can reach the theoretically optimum effi- 
ciencies no matter the optics behind the modulator, including 
the etalon. The way for calculating the required retardances 
for the two LC VRs are explained along with a number of rules 
and periodicities in the solutions. Section|4]analyzes these in- 
struments in terms of the influence of temperature and voltage 
instabilities, as well as of thickness inhomogeneities (rough- 
ness), of both the LCVRs and the etalon(s), on the final mag- 
netographic and tachographic measurements. Finally, Sect. [5] 
summarizes the results. 

2. THE THRESHOLDING ACTION OF RANDOM NOISE 

Most astrophysical measurements are nothing but photon 
counting. Their accuracy, therefore, depends on photometric 
accuracy, that is, on a battle between our ability to detect 
changes in the solar (stellar) physical quantities and the noise 
that hide such changes. The key concept is changes: we need 
to discern if a given quantity like the magnetic field strength, 
B, or the line-of-sight (LOS) velocity, vlos, varies among pix- 
els: whether or not it is greater or smaller than in the neighbor 
zones. The only tool we have to gauge these changes is 
the observable Stokes parameter changes that are linked to 
them through the response functions. Discussing response 
functions is out of the scope of this p aper as they have been 
exten sively analyzed elsewhere (e.g.. [Del Toro Iniesta et all 
2010t lOrozco Suarez & Del Toro Iniestal " 120071 

Del Toro Iniesta & Ruiz Cobol 
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As explained by iDel Toro Iniesta & Martinez Pillet 
purely phenom enological approaches (e.g., 
ICabrera Solana et"aT1 120051) are also valid to establish a 
relationship between changes in the physical quantities and 
the Stokes parameters. Here we shall concentrate on noise; 
on how it establishes the minimum threshold below which 
no signal changes can be detected. In spectropolarimetry, 
the customary estimate for noise (and thus for the detection 
threshold) is the standard deviation of the continuum signal 
because polarization is assumed to be constant at continuum 
wavelengths. Therefore, the noise is calculated either over 
a continuum window in a given spatial pixel or over all the 
spatial pixels of a map in a given continuum wavelength 
sample. Both estimates should agree as it is the case in most 
observations. If we call S = (Si, S2,Ss, S4) the (pseudo-)vector 
of Stokes parameters and denote by 07, with i = 1,2,3,4, 
the standard deviation of each Stokes parameter, then the 
signal-to-noise ratio in each parameter, (S/N)j, is defined 
as the inverse of this deviation in units of the continuum 
intensity: 



(S/N)i = 



S, 



(1) 



where index c refers to continuum. Thus, when we say 
that our observations have, for example, a S/N = 1000, we 
mean that when signals in a given (non specified) Stokes 
parameter are greater than 10~ 3 Si. c can be detected and 
this is certainly valid for that parameter but not necessar- 
ily for the others. As a matter of fact, if noise is random 
(or uncorrelated with th e signal) and can be represented by 
a Gaussian distribution dKeller & Snik 2009), according to 
IDel Toro Iniesta & Colla dos (2000|h 



where e, stands for each one of the so-called polarimetric ef- 
ficiencies of the instrument. The polarimetric efficiencies de- 
pend in a non-linear way on the modulation matrix elements 
(cf. Eq. [37l i that on their tu rn come from the first row of the 
polarimeter Mueller matrix (IDel Toro Ini esta 2003). Since all 
the efficiencies are necessarily less than the first (that of the 
intensity), Eq. (f2]i means that the noise is always larger in the 
polarization para meters than in the intensity. Then one can 
easily see that (see Martinez Pillet et al.ll 19991) 



(S/N) i =-(S/N) u i = 2,3,4, 
£1 



(3) 



that is, that the signal-to-noise ratio for Stokes 52, S3, and 
S4 is always less than that for Stokes Si. Let us point out 
here, however, that other systematic (or instrumental) errors 
like those introduced by flat fielding of images may invalidate 
the above equation. We are explicitly discarding these other 
sources of noise from our analysis. 

The Stokes parameters cannot be measured with single ex- 
posures. Instead, for vector polarimetry, a number N p > 4 
of single detector shots are recorded each providing a linear 
combination of all the four Stokes parameters. The set of N p 
individual measurements constitutes a modulation cycle that 
is characteristic of an instrument mode of operation. After 
demodulation, that is, after solving the set of N p linear equa- 
tions made up of the individual exposures, the Stokes vec- 
tor is measured. To increase the signal-to-noise ratio of the 
measurement, many instruments use accumulations, that 
is, repeat the modulation cycle N a times and the correspond- 
ing polarization images are added together. Since the degree 
of polarization of the incoming typical solar beam is fairly 
small, each single shot usually has the same light levels and, 
hence, the same (photon-noise-dominated) signal-to-noise ra- 
tio s/n. Thus, the si gnal-to-noise ratio of Si is related to the 
single-shot s/n (see iMartmez Pillel 120071 for an illustrative 
description) through 



(S/N) 1 = (s/ri)e 1 ^/N~N a , 



(4) 



Ei 

(7j = — CTi, I = 1,2,3,4, 



(2) 



because Si is often retrieved from the sum of all the accumu- 
lated images for all the polarization exposures of the given 
modulation scheme. An advisable practice for characterizing 
the signal-to-noise ratio of an instrument is to always refer 
to that of intensity and equate S/N = (S/N)\. This is conve- 
nient because there is only one intensity while the other po- 
larization Stokes parameters are three and one would need to 
specify which one is meant each time. Let us remark, how- 
ever, that this convention is not universal and some authors, 
always interested in the polarization features, think of and 
speak about some of the other three Stokes parameter signal- 
to-noise ratios. Such an alternative convention makes sense if 
one takes into account the differential character of polariza- 
tion measurements: demodulation implies that Stokes S2,3,4 
are essentially retrieved from image differences; hence, any 
systematic error like that produced by flat fielding is naturally 
mitigated (or eventually cancels out). On the contrary, the ad- 
ditive character of Stokes Si implies that intensity noise can 
be higher than simple photon noise. We shall hereafter follow 
the S/N = (S/N) 1 convention in the paper for the sake of sim- 
plicity in the description and in the equations (thus, explicitly 
neglecting systematic errors). When people is more interested 
in the other three Stokes parameter signal-to-noise ratios, Eq. 

(0 provides the obvious help. 

As demonstrated by IDel Toro Iniesta & Colladosl (120001) . 
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the maximum efficiencies that an ideal system can have are 
e = (1 , 1 /y/3, l/\/3 , 1 /V3), if all the three last ones are equal. 
Therefore, the relationship (S/N) t < l/y/3(S/N), i = 2,3,4, 
holds for any random-noise-dominated polarimetric systemQ 
Or, in terms of Eq. (fJJ, necessarily, 

a t >V3ai. (5) 

Equation (0 means that the detection threshold is bigger for 
the polarization parameters than for the intensity. Detectabil- 
ity is smaller in polarimetry than in pure photometry. 

An important parameter describing the state of any beam of 
light is its degree of polarization 



i 4 

^1 ;=9 



(6) 



Now the question naturally arises as to what is the minimum 
detectable degree of polarization by a given polarimetric sys- 
tem. If the uncertainties in the Stokes parameters are uncor- 
rected (and this should be especially true when considering a 
statistics on all the pixels of an image), error propagation in 
Eq. © gives 



(7) 



where, for convenience we have made s, ■ = Sf, i = 1,2,3,4. 
Now, since a Sj = 2Siai, it is easy to see that 




(8) 



and, according to Eq. (f2]i and since oy = 2pa p , one can write 



<? 2 F 2 \ rr 2 
3,- £j \ (7, 




(9) 



Now, if all the three last efficiencies are the same and certainly 
less than their maxima, we finally obtain 



> 



S/N 



an 



inequality 



already 



published 



(10) 



by 



(I2010I) . 

and 

without demon- 



Del Toro Iniesta & Orozco Suarez 
Martinez Pillet et all 20 1 lb. 

Del Toro Iniesta & Martinez PilleJ (120101) 
stration. 

Our instruments are aimed at measuring magnetic fields and 
velocities. Therefore, any reasonable design should include 
lower limits for these quantities within the overall error bud- 
get. Detectability thresholds for the Stokes parameters im- 
ply thresholds for the magnetograph and tachograph signals 
as well. The rest of this section is devoted to estimate them. 
Of course, any estimation that one can make depends not only 
on the instrument but also on the inference technique. Most 
modern magnetographs use inversion of the radiative transfer 

1 There may be polarimeters that are designed to measure not all four 
Stokes parameters or that aim at better accuracies for given Stokes param- 
eters. In such cases some (never all) efficiencies can be greater than 1 /v3. 
In what follows, however, we assume that our interest is the same for S2, S3, 
and 54. 



equation to infer values for both the magnetic field vector and 
the plasma velocity. These inferences involve all four Stokes 
parameters and, hence, should be more accurate than those us- 
ing just one or two of them. However, for the sake of clarity 
in the analytical derivation, we shall consider errors induced 
in the magnetographic and tachographic formulas (fill , (fl2t . 

and ([nil- 
Using classical magnetographic formulas, the longitudinal 
and transverse components of the magnetic field are given by 



B 



km ' 



1 c 

Jfc 



and 



^tran — ^tra 




(ID 



(12) 



where ki on and k tmn are (model-dependent) calibration coef- 
ficients and V s and L s are the circular and linear polarization 
signals calculated as 



1 " A 

v s = — VV&w, 



(13) 



where a, = 1 or -1 depending on whether the sample is to the 
blue (including the zero shift) or the red side of the central 
wavelength of the line, and 



1 nx 



(14) 



In the above equations, n\ stands for the number of wave- 
length samples. If we now assume that the minimum polar- 
ization signals are V s = 04 and L s = 02 = 0-3, the minimum de- 
tectable thresholds are 



S(B ]on ) = 



klon £l 

J/N el 



and 



<S(fitran) = K 



\ e\js2 
S/N 



= k„ 



l ei/ £3 
S/N ' 



(15) 



(16) 



Expressions (fl"5T l and ([Tol l give the explicit dependence of 
magnetic detectability thresholds in terms of the instrument 
efficiencies and are very useful in practice. For exam- 
ple, if we use the calibrati on constants for IMaX quoted in 
iMartmez Pillet et ail (1201 11) . assuming that the maximum po- 
larimetric efficiencies have been reached, and typical signal- 
to-noise ratios of 1700 (~ 1000 for 52,^3, and 54), the mini- 
mum longitudinal and transverse components of the magnetic 
field detectable with that magnetograph are 5 and 80 G, re- 
spectively. 

As far as the velocity is concern ed, we shall assume that 
the Fourier t achometer techn ique dBeckers & Brownl 119781 
Brown 1981; Fernandes 1992) is used: 



2c6X S\-g + S\ -3 — 5i,+3 — Sl,+9 

vlos = — — arctan - 



7rA(] 



S 



1.-9- 



>l,-3- 



>l,+3 



+ S 



(17) 



l,+9 



where c is the speed of light, <5A is the spectral resolution of 
the instrument, and Ao is the central wavelength of the line; 
-9, -3, +3, +9 stand for the sample wavelengths of the inten- 
sity, measured in picometers with respect to Ao. Let us assume 
that the minimum detectable difference between symmetric 



4 



Del Toro Iniesta & Martinez Pillet 



wavelength samples (such as Si -3 — S1.+3) due to LOS veloc- 
ity shifts is <j\ . Then, if the difference between the samples 
at the same flank of the line is approximated by 1/2 S\ yC , the 
minimum detectable LOS velocity can be approximated by 



N 2cSX 2 
(5(v L os) ^ — r— arctan— -. 

ttXq S/N 



(18) 



Likewise Eqs. ( TTOb . ( TTBT l and (fTol i. this new expression ( fTBl 
relates the velocity threshold with the signal-to-noise ratio of 
the instrument. If we use again IMaX values (SX = 8.5 pm; 
Ao = 525.02 nm) and assume S/N = 1700, the minimum de- 
tectable LOS velocity change is roughly 4 ms" 1 . 

2.1. Uncertainties induced by photon noise 

Fluctuations in the light levels due to photon statistics nec- 
essarily imply variances in the Stokes parameters that in the 
end induce uncertainties in the measured physical quantities, 
B\ n, #tran, and vlos- in this section, we are going to estab- 
lish a relationship between those variances and uncertainties. 
Note that we discard for the moment any random fluctuation 
in the instrument that will be dealt with in Sect. @] 

Error propagation in Eq. (ITU easily yields 



bL \s/n) ' 



(19) 



because a v = o~\/n\. Using now Eqs. (0, (1111) . and ([TJ, one 
obtains 



n x B 



+ 1 



Ion c 4 



1 



(S/N) 2 



(20) 



that relates the Z?i on relative error with itself, the signal-to- 
noise ratio of the observations, and the polarimetric efficien- 
cies. On its turn, error propagation in Eq. (TT2l gives 



B, 



q \s/n 



(21) 



Now, since the variances for Stokes 5? and Stokes S3 should 
be approximately the same, of ~ (1 l2n\)(o\+<ji^ and, using 
Eqs. ©, ( ITZl ). and dTJ, Eq. ( 1211 ) turns out to be 



B, 



1 



1 



(S/N) 2 



(22) 



that again relates the relevant magnetographic quantity rela- 
tive error with itself, the photon-induced signal-to-noise ratio 
of the observations, and the polarimetric efficiencies of the 
instrument. 

If we use the values for t he calibration constants quoted in 
iMartfnez Pillet et al.1 (1201 ll) for IMaX, n\ = 5 for this instru- 
ment, and assume that the maximum polarimetric efficiencies 
are reached, then the estimated relative errors for Z?i on and 
fltran induced by a photon noise of S/N = 1700 are of 2 and 
15%, respectively, for magnitudes in either quantities of 100 
G; for magnitudes of 1000 G, the relative errors drop to 0.2 
and 0.1%, respectively. 

After a similar calculation for photon-noise-induced uncer- 
tainties in the tachographic formula (flTt . one gets 

2 _ 4c 2 (SX) 2 2a\ 



n 2 X 2 



A 



(23) 



slight asymmetry between Eq. ( 1231 and Eqs. (EOt and ( 1221 is 
not such as the ratio a\ / A is a kind of inverse, square signal- 
to-noise ratio. A numerical estim ate for the IMaX instrume nt, 
and using the FTS spectrum by iBrault & Neckell (119871) to 
evaluate A for its Fe I line at 525.02 nm, we conclude that 
the photon-noise-induced uncertainty is 4 ms" 1 . 

3. AN OPTIMUM VECTOR PLUS LONGITUDINAL POLARIMETER 



As explained by lMartfnez Pillet et al.l (120041) . a versatile po- 
larimeter is obtained through the combination of two nematic 
liquid crystal variable retarders (LCVRs) with their optical 
axes properly oriented at 0°and 45°with the Stokes S2 posi- 
tive (X) direction. This is so because it can provide optimum 
modulation schemes for both the vectorial and the longitudi- 
nal (S\ ± 54) polarization analyses by simply tuning the volt- 
ages that change their retardances. The theoretical maximum 
efficiencies mentioned above can in principle be reached by 
such an ideal polarimeter. We have assumed these maximum 
efficiencies for our instruments so far. However, instrumental 
effects may corrupt the measurement so that the final efficien- 
cies are lower. Let us see in this section what happens if some 
typical optical elements are included between the modulator 
and the analyzer in the analysis. 

The corrupting effect of the optical elements of an instru- 
ment in the final polarization analysis is called instrumental 
polarization. It is well known that those optical components 
acting on light after the polarization modulation do not pro- 
duce any instrumental polarization. However, nobody has yet 
demonstrated whether optimum polarimetric efficiencies can 
still be reached no matter the optics in between the modulator 
and the analyzer. In this section we are going to show that this 
is the case with these two-LCVR-based polarimeters because 
retardances can be fine tuned by simply changing the acting 
voltages. This property certainly makes this type of polarime- 
ters very versatile and optimum for solar investigations. To 
understand the result, let us start by demonstrating that, in- 
deed, a polarimeter made up of two nematic LCVRs oriented 
as above plus a linear analyzer can reach the optimum polari- 
metric efficiencies; 

According to iDel Toro Iniesta! (120031) . the modulation ma- 
trix of any polarimetric system consists of rows that equal the 
first row of the system Mueller matrix for each of the mea- 
surements. If R(9,6) stands for the Mueller matrix of a gen- 
eral retarder whose fast axis is at an angle 9 with the X axis 
and whose retardance is 5, our LCVR Mueller matrices can 
be described by Mi = R(0,p,) and M 2 = R(7r/4,r,), where 
i = 1,2,3,4 is an index for each of the four measurements. 
Hence, in our case, where the analyzer (of Mueller matrix M4) 
is a linear polarizer at 0°Q such a mod ulation matrix, disre- 
gardi ng a 1/2 gain factor, is given by (see Martinez Pillet et al.l 
l200l 

/ 1 cosri sinpisinr; -cospisinri \ 
~_ 1 cosr 2 sinp2sinr 2 -cosp 2 sinT 2 
_ 1 COST3 sinp3sinr3 -cos/93 sinr3 ' ' 
V 1 COST4 sinp4sinr4 -cosp4sinT4 J 

As explained in IDel Toro Iniesta & Colladosl ( 120001) . if all 
the three last column elements of O have a magnitude 
of l/\/3 (with their signs properly altered), then the 
modulation is optimum and the maximum efficiencies are 



where A = (S\ _9-5i !+ 3) 2 +(5i i+ 9-5i -3) 2 , and the variances of 
the Stokes S\ samples are all assumed to be a\ . Note that the thou s h ^ double calculation is not necessary. 



2 Dual-beam polarimeters use a polarizing beam splitter as a double an- 
alyzer. Hence, another analyzer at 90° is indeed present simultaneously al- 
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reached. |cosr| = 1/V3 has the four independent solutions 
t = 54. 736,125.°264,234.°736, and 305.°264. With them, 
Isinpsinrl = l/y/3 is equivalent to |sinp| = V2/2 that has 
four independent solutions as well: p = 45°, 135°, 225°, and 
315°. The verification of the above two equations ensures 
the automatic verification of that for the third column and, 
therefore, we have found that several combination of matrix 
elements exist that qualify O as the modulation matrix of an 
optimum polarimetric scheme, as we aimed at demonstrating. 

Real polarimeters, however, have some optics in between 
the modulator and the analyzer. Very importantly, modern 
magnetographs like CRISP, VIP, IMaX, or SO/PHI have one 
or several Fabry-Perot etalons. Such etalons can modify the 
Mueller matrix that leads to a modulation matrix like that in 
Eq. (l24l and, hence, we must check whether or not the re- 
sulting modulation matrix, O', remains optimum. To do that, 
let us model the most general behavior of an etalon as a re- 
tarder M3 = R(6* et aion, ^etalon)- Then, the final Mueller matrix 
of the system is now F = M4M3M2M1 and its first row (again 
disregarding the gain factor) is given by Fu = 1, 

Fn = M 3 , 2 2 COS T, +M 3 , 24 SUIT;, 

F13 = M3 22 ship, sinr,+M3 23 cosp,— M3 24 sinp; cost,-, 
F14 = — M3 22 cos pi sin Tj +M3 23 sin p t +M3 24 cos p; cos 77 . 

(25) 

We do not need any more matrix elements of F because the 
rows of the new modulation matrix are 0[j = F\j(Ti,pi). Now, 
we only need to find out four different combinations of the 
first and second retardances that are solutions for Eqs. (|25l 
with \F\k \ = 1/Vj, where k = 2,3,4. Equations (l25T l are tran- 
scendental and, thus, have to be solved numerically. However, 
before proceeding with the numerical exercise we can realize 
several features in the solutions. First, the trivial cases, where 
<Wion = (that is, no etalon exists or it is not birefringent) 
or Etalon = 0,7r/2, the orthogonal directions of the analyzer 
axis, are indeed trivial because the effect of M3 disappears 
and O'jj = Ojj. Second, a number of periodicities can be de- 
duced from the equations structure: 



• If to is a solution for the first of Eqs. (|25> with F12 = 
l/y/3, then To + (2£+l)7r, with k integer, are solutions 
for that equation when F12 = -1/ v3 and vice versa. 



• If po is a solution for the second or the third of Eqs. 
with F\2 =1/ \/3, then po + (2k+ l)n, with k integer, are 
solutions for that equation when F12 = — 1/ V3 and vice 
versa. 



• If po is a solution for the second of Eqs. ( |25l l with 
F12 =1/ V3- then po + (2k+ V)tt/2, with k even integer, 
are solutions for the third equation when F12 = -1/ y/3. 
When k is odd, then the solution for the third equation 
is when F12 = 1 / v3 as well. 

To solve the first of Eqs. d25l l let us consider the function 

/(t) = F12 - A/3,22 cos t - A/3,24 sin r, (26) 

that has extrema where its derivative becomes zero. This oc- 
curs at to where either M3 24 = sin to and M3 22 = cos to or 
A^3,24 = - sin to and M3.22 = - cos to. These values imply that 
/(Tmax) = F12 + 1 and /(Tmin) = F\2 — 1. That is, the maximum 
of the function is positive and the minimum is negative when 
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Figure 1. Retardances of the first (bottom) and second (top) LCVRs as func- 
tions of the angle 8 of the etalon orientation with respect to the 5? positive 
direction. Only values between and 2n are displayed. Different line types 
refer to different values of the etalon retardance (see text). 

IF12I = l/v3 (which is required for reaching optimum effi- 
ciencies)0 Therefore, since / is continuous, Bolzano's theo- 
rem ensures that a solution exists in (T m i n ,T ma x) and this en- 
ables us to find that solution, for instance, through the bisec- 
tor method. This has to be done just once per value of F12; 
the other value derives from the first of the above specified 
properties. 

As a summary, the first of Eqs. (l25"t has four solutions in 
[0,27r], each two belonging to one of the signs of F12. For 
each of these four retardances, t,, four possibilities are open 
according to the values of F13 and F14. These four solutions 
for pi in the second and third of Eqs. d25l l can be shown to be 
enclosed in the following single expression, 



cos p = ± 



(n =p m) 
^(m 2 + « 2 )' 



(27) 



where n = M393 and m = M322 sin T-M324 cost. A further 
property of the solutions thus derives from Eq. (|27t : if po is a 
solution for the second and the third of Eqs. (l25l l, then lir-p^, 
ir-po, and ir + po are solutions as well. 

Therefore, the presence of an etalon modeled as a retarder 
is not a problem for the two-LCVR-based polarimeter to be 
optimum. No matter the possible retardance or orientation the 
etalon may have, we are always able to find out more than four 
combinations of p and r that ensure theoretical polarimetric 
efficiencies for all three Stokes parameters all equal to 1 / V3. 
In practice, these new solutions can be achieved by simply 
tuning the acting voltages of the two LCVRs. 

FigureQ]displays the LCVR retardances in [0,27r] that guar- 
antee optimum performance as functions of the orientation 
angle of the etalon. Values for the second LCVR are in 
the top panel and those for the first one are in the bottom 
panelQ Different colors correspond to different solutions; 
different line types correspond to different values of the as- 
sumed etalon retardance: 0°(dotted), 30°(solid), 45°(dashed), 
and 60°(dashed-dotted). As commented on before, when 

3 Note that we have selected one out of the infinite solutions for the deriva- 
tive of / to be zero, but this is coherent with our neglecting multiplicative, 
gain factors in the definition of Mueller matrices. 

4 Note that retardances larger than 27r may be needed for design conve- 
nience. Their values are easily deducible from the above-mentioned proper- 
ties. 
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#etaion = 90°, both p and r recover the same value as if the 
etalon were absent. Moreover, both retardances are periodic 
with 6, with periods tt/A (p) and ir/2 (r). It is also interesting 
to note that four out of the eight solutions for p are equal to 
the other four but phase shifted by n/4. 

Now, it is a little tedious but easy to demonstrate that, re- 
gardless of how many, mirrors can be introduced in the optical 
path between the modulator and the analyzer (as in real instru- 
ments) without affecting the maximum polarimetric efficien- 
cies, provided they all are perpendicular to the optical axis 
plane. From the Mu eller matrix of a single mirror, one can 
realize (ICollettl fl 992h that the matrix of such a mirror train, 
no matter the angles between them, keeps always the shape 



E = 




(28) 



The elements of the first row in the new final Mueller ma- 
trix of the system become F[ - = {a+b)F\j, j = 1,2,3,4. That 
is, the final modulation matrix remains the same as before 
introducing the mirrors but scaled by a gain factor that can 
be disregarded as we have been doing for all the treatment. 
Therefore, we can conclude that optimum efficiencies can still 
be achieved with as many mirrors as needed. Since a mirror 
is indeed a combination of a retarder and a p artial polarize r 
(the d element is zero for the latter; see e.g. iStenflol [19941 
the same conclusion can be reached for whatever differen- 
tial absorption effects for the orthogonal polarization states 
that may be located between the modulator and the analyzer. 
Therefore, if, for example, the etalon or the diffraction grating 
of the instrument display different transmittances for orthog- 
onally polarized beams the polarimetric efficiencies can still 
attain maximum values. 

4. INSTRUMENT-INDUCED INACCURACIES 

Now that we know that our magnetograph can reach opti- 
mum polarimetric performance, let us study the behavior of 
this particular instrument against instabilities in its main opti- 
cal elements. Photon noise is not the only harm for magneto- 
graphic or tachographic measurements. Instabilities of differ- 
ent types like those in the temperature or in the tuning volt- 
age of both the LCVRs and the etalon, or roughness in their 
final thicknesses, can induce inaccuracies. For single mea- 
surements, the inaccuracies can imply errors in magnetic field 
or absolute wavelength calibration. For time series like those 
needed in helioseismological studies, such inaccuracies may 
avoid detection of some particular oscillatory modes. An as- 
sessment on such inaccuracies is therefore in order for clearly 
defining design tolerances of these instrumental quantities. 
Such tolerances should ensure the fulfillment of the scientific 
requirements of the instrument. 

4.1. Polarimetric inaccuracies 

An important example of a scientific requirement is the po- 
larimetric accuracy of the system. By such we understand any 
of the inverse signal-to-noise ratios for S2, S3, or S4 as defined 
in Eq. (HJ. Following our general assumption that these (S/N), 
are intended to be the same by design, aiming at a S/N of 
1700 is equivalent to require the system to have a polarimet- 
ric accuracy of 10~ 3 . Temperature, voltage, and other insta- 
bilities and defects of the LCVRs lead to changes in the retar- 
dances that, on their turn, induce modulation and demodula- 
tion changes. Such changes can be seen as cross-talk between 



the Stokes parameters that drive to covariances in the magne- 
tographic measurements (lAsensio Ramos & Colladosl 120081 : 
se e also an interestin g discussion on seeing-induced cross-talk 
in lCasini et alJl201 lb . Since our modulation matrix (l24l is an- 
alytical, we attempt an analytical approach to the study of the 
effect of these retardance changes onto the polarimetric accu- 
racy of the system. 

Let D be the demodulation matrix of the instrument that 
always exists. Thus, OD = DO = 1. The measured Stokes 
vector is then given by S = DI meas , where I meas stands for the 
four intensity measurement vector of each modulation cycle. 
If we linearly perturb matrix D as a consequence of a small 
perturbation in the LCVR retardances (<r Pt ,a Tk , k= 1,2,3,4), 
then the Stokes vector becomes S + S' = (D+D^I mea s, where 
the perturbed demodulation matrix D' is given by 



1V W^D 9D 



(29) 



It is obvious that the polarimetric accuracy requirement di- 
rectly implies that none of the four elements of S' = D'I meas 
can be greater than 10~ 3 . To fulfill that requirement, let us 
then study how the perturbation in the retardances can be pro- 
duced and which are the tolerances for the instrument quanti- 
ties whose fluctuations produce them. 

If we call §l anyone of the retardances, we have by defini- 
tion that 



An 



(30) 



where f3 = n e — n is the birefringence, i.e., the difference be- 
tween the extraordinary and the ordinary refractive indice s of 
the liquid crystal, and t stands for its geometrical thickness. 
Therefore, it is evident that 



a 



5 L 



S 2 L /3 2 f 2 



(31) 



that is, the relative inaccuracy in the LCVR retardance is the 
square root of the sum of the square relative inaccuracies in 
the birefringence and in the geometrical thickness at the oper- 
ating wavelength. Let us ascribe for convenience any possible 
local fabrication defect in the LC like an air bubble to the 
thickness inaccuracy, so that birefringence can be considered 
spatially constant for the whole device. 

Variations in the birefringence can be produced by either 
variations in the LCVR temperature, the acting voltage, or 
both. Hence, one can write 



2 22 22 
Og =q T a T + q v a Vl 



(32) 



where qj and qy are values of the (partial) derivatives of j3 
with respect to T and V at the given values of voltage and 
temperature, respectively. Since we indeed have calibrations 
of the Sl dependences rather than those of f3, we better rewrite 
Eq. (ED as 



X2 

2 2 2, 2 2 , °L 2 

<t Sl =m T a T + m v a v + -^a t , 



(33) 



where m\ and m v have a clear meaning and can be deduce d 
from calibration s. According to [Martinez Pillet et ail (1201 11) . 
based on data bv lHeredero et al.l (12007b . 



m r =-1.16 + 0.305V-0.02V 2 



(34) 
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Figure 2. Retardance (in degree) of a specific nematic LCVR as a function 
of the acting voltage (in volts; solid line) and its derivative with inverted sign 
(dashed line). 

for V < 8 volt and otherwise, and m v can be obtained from 
data like those displayed in Fig.|2]where 8l (solid line) and its 
derivative (dashed line; inverted sign) are plotted as functions 
of the acting voltage for a particular LCVR. 

To get a numerical estimation of the value of real tolerances, 
that is, of the maximum ay, cry, and a, affordable in real in- 
struments in order not to have the S'-elements greater than 
10~ 3 , we shall use IMaX parameters^plus the analytic expres- 
sions for D and its partial derivatives^ The last term in Eq. 
d33l may vary spatially and is, thus, responsible for the pixel- 
to-pixel variations of the retardance but can easily be cali- 
brated if needed. Indeed, roughness in the device thickness is 
a fabrication specification and can be checked upon delivery 
from the manufacturer. Thickn ess inaccuracies may produc e 
locally significant effects (e.g. jAlvarez-Herrero et al.ll2010t) . 
Our own estimation, using Eqs. d29l [33l indicate that relative 
errors in the thickness larger than 6 % induce perturbations of 
the Stokes vector that are larger than the polarimetric accu- 
racy. Therefore, should this be the only instrumental instabil- 
ity, specific pixel-to-pixel calibration of D would be needed if 
the relative roughness is larger than 6 %. Note that, since the 
typical thicknesses of LCVRs are of the order of micrometers, 
a roughness less than 6% may mean a stringent requirement 
for the manufacturer of the order of tens of nanometers. 

Using again Eqs. (l29l 1331) . we find out that instabilities 
larger than oj = 600 mK or ay = 1 .5 mV deteriorate the po- 
larimetric accuracy below the required 10~ 3 . These two tol- 
erances and that for the roughness have been calculated after 
assuming that each instability is acting individually. Accord- 
ing to Eq. Q31 l. the final uncertainty in the retardance stems 
from the three sources simultaneously. Hence, a safety reduc- 
tion factor of a/3 (assuming all the three contribute the same) 
is advisable. Therefore, in the end, the final tolerance specifi- 
cation for our instrument to reach a polarimetric accuracy of 
10~ 3 is 300 mK for temperature, 1 mV for voltage, and a 4 % 
for LCVR roughness. 

5 The combination for IMaX retardances was p = [315,315,225,225] 
and r = [305.264,54.736,125.264,234.736], in degrees. Their corre- 
sponding voltages were V(p) = [2.535,2.535,3.112,3.112] and V(t) = 
[2.4,9.0,4.3,2.9], in volts. 

6 An IDL program that evaluates such analytic expressions is available 
upon request. Direct (by hand) evaluation is so tedious that only with the help 
of software applications like Mathematica such expressions can be obtained. 



4.2. Magnetographic inaccuracies 

The retardance perturbations of Eq. (03) induce changes 
in the maximum polarimetric efficiencies of the instrument. 
Such changes necessarily imply modifications in the magne- 
tographic measurements. The modifications might jeopardize 
the quality of the results. Imagine, for instance, that a require- 
ment on the repeatability of B\ on and Btran applies because 
we are interested on a measurement time series: a calcula- 
tion of tolerances in the instrument parameters (T, V, rough- 
ness, etc.) that ensure the fulfillment of the magnetographic 
repeatability is in order. 

Since no explicit dependence of Bi on and Z? tran on e, exists, 
we cannot analytically gauge the induced inaccuracies in any 
magnetographic measurement. Nevertheless, we can obtain a 
hint on the global by studying the specific variations of S(B\ on ) 
and <5(B t ranX the minimum detectable values of such magneto- 
graph quantities. 

Error propagation in Eqs. (fT31 and (ITol readily gives 



2 

a S(B lm ) ' 



S 2 (B ]on ) 



and 



a MB lmn ) : 



8\B 



\6ei, 



ran) 2 
- C 2 



5\B lmn ) 
\6e\ 



(35) 



(36) 



According to iDel Toro Iniesta & Colladosl (120001) . the maxi- 
mum polarimetric efficiencies that can be reached by any sys- 
tem are 



O 2 



(37) 



where Oj, are the matrix elements of O t , the transpose of O. 
For a system with a modulation matrix like that in Eq. (l24l . 
it is easy to see that the efficiency inaccuracies ensuing the 
thermal instabilities are 



= 0, 



(38) 



as a consequence of using normalized Mueller matrices, and 



= i^E sin22r >4' 



(39) 



7=1 



1 f 4 

2 1 ) \ ~* • 2 ^ -4 2 

^L,., = i6 I 2^ sin 2 PJ sin T J a ' Pi + 
J2 sil)4 Pj sin 2 27-^ cT^ | , 

CT L 4 = ^|E sin22 ^ sin4 ^< + 



cos 4 pj sin 2 2tj a 2 .. 



(40) 



(41) 



where err and a\ are the variances of the LCVR retardances 

~j Pi 

and where we have neglected the influence of possible time- 
exposure differences or instabilities like those caused by a 
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rolling-shutter detector or a non-ideal repeatability of a me- 
chanical shutter. Such influences can be estimated separately 
for the specific instrument and directly scale the effective ex- 
posure time (either per pixel or per frame). 

Using the values for IMaX, and assuming maximum effi- 
ciencies, instabilities of 0.3 K or of 1.1 mV produce a 5 % 
repeatability error in the threshold for Z?i on . A comparison of 
Eqs. (l35l l and (l36l l readily tells that the effect is 2 times smaller 
on the relative repeatability error in the threshold for B tlan but 
since the threshold itself is 16 times larger, the absolute error 
is in the end 8 times larger as well. 

4.3. Velocity inaccuracies 
Error propagation in Eq. (TTTt yields 



2 / dvLOi 



2 9 / dv LOS s 



dv 



LOS 



OS 



11 



2 

°1,* 



(42) 



with k = -9 , -3 , +3 , +9 . Uncertainties in the spectral resolution 
come from uncertainties in the etalon spacing and related fab- 
rication details that produce etalon roughness. Errors in the 
central wavelength come from the etalon tuning that mostly 
depends on the ambient temperature, T, and on the tuning 
voltage, V: 

ol =kW T + kWv, (43) 

where kj and ky are constants that give the (linear) depen- 
dence of Ao on T and V. Finally, uncertainties in the Stokes 
S\ samples come both from pure photon noise, as in any pho- 
tometric measurement, and from etalon tuning uncertainties. 
Since the (inexplicit) dependence of Sj^ on Aq is non linear, 
let us linearize it (that is, introduce a small perturbation and 
take the first approximation) and write 



2,2 2 



(44) 



where we have assumed that the photometric contribution is 
equal for all the samples and indeed equal to the photon noise 
as calculated in the continuum; sif stand for the derivatives 
of the Stokes S\ profile at the corresponding wavelengths. 

After a tedious but straightforward algebra, Eq. d42l can be 
recast as 

"LOS 2 , 4cW 2<7? 



^LOS 



(S\) 2 



A 



^ os 4c 2 (S\f d,+d 2 
\ 2 + 



r 2 M 



A 2 



(kjaj + kyCFy), 



(45) 



where A is defined in Sect. 12.11 d\ = (Si l+ 9 — Si -3) (sf _ 9 + 
s 2 +3 ) and d2 = (S\ i _9-5'i i+ 3) 2 (S[ +9 +s 2 _ 3 ). Hence, clear contri- 
butions to the final LOS velocity uncertainty can be discerned 
from the etalon roughness, the photon noise, the etalon tem- 
perature instability, and the etalon voltage instability. 

Quantitative estimates of the various te rms in Eq. (|45) 
can be made by using the FTS spectrum by iBrault & Neckell 
dl987l) to evaluate A, di , and c/2 for a given spectral line and 
a given instrument. Let assume the HMI and SO/PHI Fe I 
line at Ao = 617.3 nm and a spectral resolution of the etalon 
of 5X = 10 pm. The first term has a clear impact on the tacho- 
graphic results: the etalon relative roughness is directly trans- 
lated into the same vlos relative uncertainty. In other words, 
we cannot expect better accuracy in the line-of-sight velocity 
(when measured with the Fourier tachometer formula) than 



that limited by the etalon relative roughness; this means that 
a mere 0.1 pm, rms resolution uncertainty induces 10 ms" 1 
errors for speeds of 1 kms" 1 . The second term coincides with 
the righ t-han d side of Eq. (1231 and has been discussed already 
in Sect. 12.11 Since the ratio between the second and the first 
terms within brackets is of the order of 4-10 9 for velocities 
of up to 5 kms" 1 , it is the second one what really matters in 
the estimation; this means that the dependence of the profile 
shape on the central wavelength of the line is really impor- 
tant. If we use the same IMaX values of kj = 2.52 pm/K and 
k v = 3.35 10" 2 pm/V for the SO/PHI etalon, Eq. g2> gives in- 
stabilities of the order of 3.3 mK or 0.25 V induce the same 
LOS velocity uncertainty of 4 kms -1 quoted above for pure 
photon noise. Another way of seeing the same effect can 
be explained by saying that a 100 ms -1 uncertainty is pro- 
duced by either a 45 mK or a 3.4 V instabilities. These uncer- 
tainties are really important when stability during given peri- 
ods of time of the instrument is required as for helioseismic 
measurements. For single shots, uncertainties in temperature 
or voltage imply thresholds for accurate absolute wavelength 
(velocity) calibration. The importance of having included the 
measurement technique in this error budget analysis is clear: 
should one have simply used the ky and kj calibration con- 
stants above and the Doppler formula an uncertainty of just 
55 ms" 1 would have been obtained. Hence, the uncertainty 
would have been underestimated by a factor almost 2. 

5. SUMMARY AND CONCLUSIONS 

An assessment study on the salient features and properties 
of solar magnetographs has been presented. An error budget 
procedure has been followed. Special care has been devoted 
in including photon-induced and instrument-induced noise as 
well as specific measurement technique contributions to the 
final variances. We have first discussed the effect of ran- 
dom noise in the measurements and deduced useful formu- 
lae -general for every device- that provide some minimum 
detectable parameters like the degree of polarization of light, 
the longitudinal and transverse components of the magnetic 
field, and the line-of-sight velocity. The detection thresholds 
are given as functions of the polarimetric efficiencies of the in- 
strument and of the signal-to-noise ratio of the observations. 
(As a proposal, we have suggested as well to use the S/N 
for the Stokes intensity as the signal-to-noise ratio for the 
instrument.) When the random noise is photon-induced, we 
have calculated as well the relative uncertainty in the magne- 
tographic and tachographic quantities. Secondly, an analysis 
is presented for those instruments based on two nematic liq- 
uid crystal variable retarders as a polarization modulator and 
a Fabry-Perot etalon as the spectrum analyzer. Although spe- 
cific for these magnetographs, the methodology can easily be 
followed by others in order to characterize their capabilities 
and accuracies. We have demonstrated that this type of in- 
strument can indeed reach theoretical maximum polarimetric 
efficiencies because solutions always exist for the retardances 
of the two LCVRs that ensure such efficiencies, hence op- 
timizing the detection thresholds and the relative uncertain- 
ties. Very remarkably, the existence of such solutions is in- 
dependent of the optics that is in between the polarization 
modulator and the analyzer. Neither retarders nor partial po- 
larizers or mirrors (the most commonly used devices) alter 
that property. The LCVR optimum retardances do depend in 
such pass-through optics but can be fine-tuned according to 
the polarizing properties of the optics. A number of rules and 
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periodicity properties of the required retardances have also 
been deduced. These polarimeters have modulation and de- 
modulation matrices that are explicitly calculated through an 
IDL procedure that is available upon request. Thirdly, error 
propagation has yielded equations relating the variances of 
the measured Stokes vector and the solar physical quantities 
and instrument parameters, hence providing a bridge between 
scientific requirements and instrument design specifications. 
The analytic character of this particular type of instrument has 
also allowed the quantitative estimation of the mentioned un- 
certainties. Hopefully, the discussion presented in this paper 
excites (and helps) further diagnostics of other instruments. 
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AYA2009-14105-C06 and AYA201 1-29833-C06, and Junta 
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percentage from European FEDER funds. 
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